Absence of epidemic thresholds in a growing adaptive network 
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The structure of social contact networks strongly influences the dynamics of epi- 
demic diseases. In particular the scale- free structure of real- world social networks 
allows unlikely diseases with low infection rates to spread and become endemic. 
However, in particular for potentially fatal diseases, also the impact of the disease 
on the social structure cannot be neglected, leading to a complex interplay. Here, we 
consider the growth of a network by preferential attachment from which nodes are 
simultaneously removed due to an SIR epidemic. We show that increased infectious- 
ness increases the prevalence of the disease and simultaneously causes a transition 
from scale-free to exponential topology. Although a transition to a degree distribu- 
tion with finite variance takes place, the network still exhibits no epidemic threshold 
in the thermodynamic limit. We illustrate these results using agent-based simula- 
tions and analytically tractable approximation schemes. 
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I. INTRODUCTION 



Throughout human history epidemic diseases have been a constant threat. The plague of 
Athens wiped out 25-35 percent of the city's population as early as 430 BC jl|. A bubonic 
plague epidemic killed between 75000 and 100000 inhabitants of London between April 1665 

n 

and January 1666, as its population was about 460000 [1[. Smallpox became a constant 
threat to Europe throughout the eighteenth century, where growing cities like London were 



especially vulnerab 
hosts to the virus 



e to the disease on account of providing incessant susceptible immigrant 
l|. After a brief respite during the 20th century epidemics are on the 



rise again. For instance in 2010, around 655000 people died of malaria 2}]. Combating such 
epidemic diseases efficiently in the mega-cities of the future is likely to require a heightened 
understanding of the dynamics of the diseases and the social networks on which they spread. 

The horrendous examples above illustrate that epidemic diseases can have devastating 
effects on populations. It is easily conceivable that this has a drastic effect on population 
structure and specifically on the structure of social contact networks in the population. 
Simultaneously, the fate of epidemics strongly depends on properties of the contact networks, 
leading to a complex interplay. 

In the past two decades, epidemic spreading has been extensively studied on different 
complex networks to understand the influence of the social contact network structure on 

fin 

the disease prevalence [3H5|- In particular, it has been shown that a crucial determinant 
of epidemic spreading is the degree distribution of the contact network, i.e. the probability 
distribution of finding an agent with a given number of contacts. If the variance of this 
distribution is finite then there is generally a threshold of infectiousness, which the disease 
has to exceed to persist in the network |6j-l8[. By contrast, if the variance of the degree 
distribution is infinite, such as in certain scale-free networks, then the epidemic threshold 
vanishes [7Hll| and diseases can persist by hierarchical spreading from hub nodes of high 
degree [12|. Thus scale-free networks allow unlikely diseases with low infectiousness to spread 
and become endemic. 

Today, the structure of social networks is recognized as a key factor with direct impli- 
cations for epidemic dynamics and potential counter measures 13l-ll7j]. This insight has 



40 motivated the integration of real world network data into epidemic models 18M22] . Further- 



more, theoretical models have been extended by including several properties of real world 
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hierarchy 



2n 



networks such as degree constraints 231 ]. clustering 24 H26L in formation filtering 27J, social 



and nonuniform transmission probabilities 
A relatively recent addition is to consider also the feedback of the epidemic on the social 
network structure, which can be indirect, e.g. by triggering behavioural changes of agents 
0, Sj; or direct by removing agents due to hospitalization, quarantine or death. 

Modelling the network response to an ongoing epidemic leads to an adaptive network, a 
system in which a dynamical inter play between the dynamics of the network and the dynam- 
ics on the network takes place 3l|, |32J. Epidemics on adaptive networks can exhibit complex 



emergent dynamics 30 



33- 



e.g. sustained oscillations, bistability, and hysteresis) and 



emergent topological properties [30j, |33|, |37J (e.g. heterogeneous degree distributions and 
assortative degree correlations). Furthermore, the study of the social response to epidemics 
is interesting from an applied point of view because it could enable enhanced vaccine control 
and effective quarantine strategies [391 ] . 
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371 ] , direct topological feedback 



In comparison to social responses to the epidemics |30j, 
via the removal of nodes has received less attention. Previous works considered the case 
where network growth and death processes are balanced and thus the population stays 



in equilibrium and fluctuates around a fixed system size [40|, |4l|. The case of continuous 
growth was studied in 45[, where new nodes attach preferentially to high degree non- infected 
nodes. It was observed that a transition from a scale-free topology to an exponential one 
takes place as the infectiousness decreases. In another study, network growth and node 
removals have been incorporated in a single model 42|. However, the authors focused on 
the epidemic oscillations and did not consider topological effects in detail. Another related 
work focused on the interplay between network growth and dynamical behaviour in the 



431 ] . where new players preferentially attach to those 



context of evolutionary game theory 
receiving higher payoffs. 

Here, we consider the growth of a network by preferential attachment from which nodes 
are simultaneously removed due to a susceptible-infected-removed (SIR) epidemic. The 
appeal of this model lies in its paradoxical nature, in absence of the disease preferential 
attachment will lead to the formation of scale-free topologies in which the epidemic threshold 
vanishes, such that the disease can invade. However, an established SIR disease will quickly 
infect and remove nodes of high degree such that the variance of the degree distribution 
is decreased and epidemic thresholds reappear, potentially leading to the extinction of the 
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disease. 

75 One can attempt to guess the outcome by the following (slightly naive) line of reasoning: 
Observing a scale-free topology implies that the epidemic is extinct. But an extinct epidemic 
implies scale-free structure and hence vanishing epidemic threshold precluding extinction. 
Logically, the only possible solution is that the epidemic persists (unconditionally) in a 
network that is not scale free. In other words, one would expect that the coevolution of 

so epidemic state and network structure should lead to a vanishing epidemic threshold in an 
exponential network. 

One could argue that the naive reasoning is wrong because the paradox above can be 
resolved temporally, such that the epidemic goes extinct while the network is exponential, 
whereas a disease-free scale-free network develops at later times. Subsequently the network 
85 cannot be re-invaded by the disease because all infected have been removed. However, this 
temporal resolution is only feasible in finite networks. In the thermodynamic limit it can 
be easily shown that a finite number of infected survive even below the epidemic threshold, 
which precludes complete extinction and thus is sufficient to reignite the epidemic at a later 
stage. 

go In the remainder of this paper we present a detailed dynamical analysis of the epidemic 
model showing that the naive reasoning presented above is actually correct in the thermo- 
dynamic limit. Here, high infectiousness elevates the prevalence of the disease but at the 
same time limits epidemic spreading by causing a transition from scale-free to exponential 
topologies. Thus a balance is reached where the finite variance of the degree distribution 

95 is matched to the finite infectiousness, such that the epidemic persists. We use a novel 
analytical approximation scheme to understand this balance in detail, and also identify a 
small parameter region, where, after all, a temporal resolution of the paradox is observed. 



II. MODEL 



We study the spreading of a susceptible-infected-removed (SIR) disease [44| on an evolving 
ioo network. In this network a given node is either susceptible (state S) or already infected with 
the disease (state I). We start with a fully connected network of m nodes and consider three 
dynamical processes: a) the arrival of nodes, b) disease transmission, and c) the removal of 
nodes. 
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In the following we measure all rates per capita, including the arrival rate. This implies 
105 that larger populations have a proportionally larger influx of agents. This assumption is 

necessary to keep the model well-defined in the thermodynamic limit and appears plausible 

e.g. for growing cities, where the attractivity increases with size. 

New nodes arrive in the population at a constant rate q and are already infected with the 

disease with probability w. Arriving nodes immediately establish links with m of the nodes 
no selected according to the preferential attachment rule [46] : A new node establishes a link 

with an existing node of degree k with probability kPk/{k), where P k denotes the degree 

distribution (the probability that a randomly picked node has degree k) and (k) = ^2kP k 

is the mean degree. 

Disease transmission occurs clt Tcttc p on every link connecting a susceptible and an 
us infected node. Therefore, nodes with higher degree are proportionally more likely to catch 

and spread the disease. 

Removal of infected nodes take place at rate r. Because we describe a fatal disease from 

which recovery is not possible, removed nodes and their links are entirely deleted from the 

simulation and do not re-enter at a later stage. 
120 Unless mentioned otherwise, the following set of parameter values is used throughout the 

paper: m = 6, m = 5, q = 0.01. In agent-based simulations the network is simulated until 

N reaches 10 7 (the scale of present mega cities), or the network has stopped growing and 

the time reaches 10 4 . 



III. ANALYTICAL TREATMENT 



125 The dynamics on and of complex networks can be captured by a set o 
differential equations, in so-called moment-c 
develop a heterogeneous node approximation 



osure approximations 



37 



coupled ordinary 



47 



531 ] . Here, we 



7H9J, where the network evolution is captured 
in a set of equations for the node densities in different degree-classes. We define [A k ] as the 
density of nodes of state A e S,I and degree k among all nodes in the network. Furthermore, 
130 the total density of S-nodes is denoted as [S], [S] = ^2 k [Sk\- The density [I] is defined 
analogously such that [S] + [I] = 1. We note that the number of nodes N is not conserved. 
Therefore, node densities are renormalized due to node arrival and removal events. 

We now formulate a dynamical system that captures the time evolution of the densities 
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[Sk] and [Ik] due to infection, node arrival, and node removal processes, 

^[S k ] = q(^(l - w)S k>m + k[S k ) + (k- l)[S*_i]) - [S k ]^j -p Zl [I]k[S k 



^r[Ik] = q(wS k>m +^r(-k[I k ] + {k-l)[I k - 1 ]) -[I k ] \ j-pzi [I]k[S k ] 



+rU[i]((* + l)[ 

dt- j q { w5k ' m + Wj( 

+r (z T [I] ((* + l)[/ fe+1 ] - k[I k ]j + [I][h] - [I k ]\ , < k < k max (1) 

135 For understanding the equation governing the evolution of the density [S k ], consider that 
new nodes with degree m arrive at the rate q and have state S with probability 1 — w. 
Thus, the density [S m ] increases at the rate q(l — w). A newly arriving node builds a link 
to a node in the S k class with the probability k[S k ]/ (k) and causes it to pass into the S k +i 
class. Because m such links are established by each newly arriving node, the density [S k ] 

uo decreases by qmk[S k ]/(k). Similarly, nodes in the S k -i class pass into the S k class at the 
rate qm(k — l)[S k -i]/ (k) . Additionally, we need to renormalize the density [S k ] when a node 
arrives. This corresponds to a loss of the [S k ] density of q[S k ]. 

At rate p, nodes within the S k class become infected through their links with infected 
nodes causing them to pass into the I k class. The density of such links can be approximated 

us to good precision by Z/[/]A;[iSfc], where Zi = (hi)/ (k). 

Finally, nodes within the S k class pass into the S k -i class due to the removal of their 
infected neighbours. Given the density of such links, z/f/J/cfSy, the density of [S k ] decreases 
by rzi[I]k[Sk]- Similarly, nodes in the S k +i class pass into the S k class corresponding to a 
gain of rzi[I](k + l)[(Sfc + i]. As infected nodes are removed, the density of all degree classes 

150 increases due to the renormalization leading to a gain of r[I][5fc] for the density [Sk]- The 
rate equation for the density [I k ] is constructed analogously. 

In the following we refer to equation (JT]) as the heterogeneous approximation. The main 
drawback of such heterogeneous approximations is the high dimensionality of the system of 
equations, which complicates the analytical solution and thus typically necessitates extensive 

155 numerical studies. Furthermore, since it is not possible to numerically integrate an infinite 
dimensional system of differential equations, we need to introduce a degree cut-off k max by 
assuming YlT=k max +i ^ k ^ Sfc=cT ^ ne higher the degree cut-off, k max , the more precise 
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the heterogeneous approximation becomes. 

In the following, we develop a low dimensional approximation by summing over the degree 
wo classes. We consider the susceptible proportion of the population [S], the mean degree (k), 
and the mean degree of susceptibles (ks) which evolve according to 

k 

£<*> = £*^(P*] + fr]), 

k 

Tt {ks)[s] = ^i k[Sk] - (2) 

k 

Using equation ([1]) and equation ([2]), we obtain 
±[S) = g (i - w - [S]) - p { -^^[S][I] + r[S)[I), 

±(k) = q(2m - (k)) + r(2(k s )[S] - (k)(l + [S])), 

d , h \ _ (k s )) (k s )\ (ki)[I} (( 2) , , 2 . (fcg)^/) 

In the following, we refer to equation fl3]) as the coarse-grained heterogeneous approxima- 
tion. 

165 Because we have not derived an equation for the second moment of the susceptible degree 
distribution (kg), equation (jHJ) does not constitute a closed dynamical system. We address 
this problem by replacing (k s ) by (k$) 2 + (ks) in an additional approximation. We note that 
this approximation is valid exactly when the network has a Poissonian degree distribution. 
It can therefore be thought of as a 'random-graph approximation'. This approximation will 

170 certainly fail in the case of scale-free degree distribution because of the degree distributions 
with diverging variance, i.e. (k s ) — > oo. By contrast, as will become apparent below, the 
approximation still performs well for distributions with large finite variance. 
Using the random-graph approximation we obtain 

i[S] =q{l-w- [S]) -P { -^I^[S}[I] + r[S}[I], 
= q(2m - (k)) + r(2(k s )[S] - (k)(l + [S])), 

d = q ((i- w )( m -(k s) ) + jks) \ _ p (k s )(k m r (MM [Ih (4) 



at V [S] (k) J (k) (k) 

where [J] and (kj) are given by the conservation laws, such that the system constitutes a 

175 closed model. In the following we refer to this model as the homogeneous approximation. 
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IV. BASIC ANALYSIS 

In this section, we present a detailed analysis of the analytical model and confirm the 
results by comparison with agent-based simulations of the network. 

Before we launch into a detailed discussion of the model, let us consider the limiting case 
iso of network evolution in the absence of the epidemic. In this case the model is identical 



to the Barabasi-Albert model of network growth [461 ]. which is known to lead to scale- free 
topologies, where the degree distribution follows a power law oc k^ 1 with exponent 7 = 3. 
In the present model the emergence of scale- free topologies is thus expected in the limit where 
the disease goes extinct or remains limited to a finite number of infected nodes <C N. When 
185 the epidemic is present, high degree nodes are disproportionately likely to become infected 
and subsequently removed, which can be expected to prevent the formation of scale-free 
topologies. 

We confirm this intuition by plotting degree distributions for various parameter sets in 
figure [U Because the homogeneous approximation does not provide any information on the 

wo degree distribution, we show a comparison of the heterogeneous approximation with agent- 
based simulations. The figure shows a good agreement between the modelling approaches 
and confirms basic intuition. When all arriving nodes are susceptible (w = 0), a scale- 
free degree distribution with the expected exponent 7 = 3 is formed for p = 0. At finite 
infectiousness p, the topology changes from scale-free to exponential. The same behaviour 

195 is observed at higher rates of infected arrivals, < w < 1. 

When all arriving nodes are already infected (w = 1), the distribution has a bimodal form 
for p = with high degree contribution coming from the initial susceptibles which never get 
infected. At positive infectiousness p, these individuals eventually die and the mode at high 
degrees disappears. 

200 In order to quantify the topological transition from the scale-free to the exponential degree 
distribution, we plot the variance a 2 of the degree distribution as a function of infectiousness 
p and fraction of infected arrivals w in figure EJ As either parameter increases, the disease 
prevalence in the steady-state, [/]*, increases and removal occurs at a high rate. As a result 
the degree distribution becomes narrower and the degree variance decreases. It is apparent 

205 that very high values of the variance are only found for low infectiousness p, whereas higher 
infectiousness quickly lead to narrow distributions. 
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FIG. 1. Degree distribution for varying infectiousness p and infected arrival fraction w. Left: 
individual-based simulations. Right: heterogeneous approximation. Scale-free degree distributions 
are observed when disease infectiousness vanishes (p = 0) . For increasing infectiousness the degree 
distribution quickly becomes exponential. Parameters: r = q = 0.01, tuq = 6, m = 5. Shown are 
averages over 1000 simulation runs. 



Above we computed the variance a 2 of the degree distribution of networks. One concern in 
any computation of this kind is finite-size effects. In the heterogeneous approximation these 
effects appear directly in form of the maximal degree that is considered in the approximation. 
210 In agent-based simulation a similar cut-off exists as the maximal degree in a network of 
finite size is bounded by the number of nodes. Hence all moments of the degree distribution, 
including the variance, must be finite regardless of the shape of the degree distribution. 
However, if the finite networks are drawn from an ensemble that becomes scale-free in the 
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FIG. 2. Variance a 2 of the degree distribution. Degree variance a 2 is plotted as a function of 
the infectiousness p and the infected arrival fraction w from agent-based simulations. The degree 
variance a 2 decreases with increasing infectiousness p and fraction of infected arrivals w. Inset: 
Degree variance a 2 as a function of the degree cut-off k c . A transition from a cut-off dependent 
(p = 0, scale-free) to an independent (p = 0.0004, p = 0.001, and p = 0.005, exponential) regime is 
observed as the infectiousness p is increased. In network simulations (circles), nodes were restricted 
to at most k c neighbours. In the heterogeneous approximation (solid lines) the cut-off k c is directly 
imposed as k nmx . Parameters: r = q = 0.01, tuq = 6, m = 5, k c = 5 x 10 3 . 

thermodynamic limit the variance a 2 is often found to increase logarithmically with the 

2i5 imposed cut-off [23 ] . 

We now rule out that low values of the degree variance a 2 , observed above, were due 
to finite-size effects by considering the variance a 2 as a function of the degree cut-off k c 
(figure [2] inset). For the case of p = 0, where we observed scale-free behaviour, we find that 
the observed variance increases logarithmically as expected. Conversely, for the finite values 

220 of infectiousness p the observed a 2 is insensitive to a sufficiently large cut-off. In summary 
these results show that fatal diseases should relatively quickly destroy scale-free structure of 
social networks at all but the smallest removal rate and/or infectiousness. 

Let us now investigate the effect of the emergent network topology on the prevalence of the 
disease. Plots of the disease prevalence as a function of the infectiousness p, and the fraction 
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FIG. 3. Dependence of the disease prevalence [I]* on infectiousness p (top) and fraction of in- 
fected arrivals w (bottom). The disease prevalence increases monotonically with infectiousness 
p and fraction of infected arrivals w. In agent-based simulations (circles), [I]* is calculated over 
the surviving runs among 10 3 total realizations. Homogeneous approximation (solid lines) is the 
analytical solution of equation Q. Heterogeneous approximation (dashed lines) is the stationary 
value of the numerical integration of equation ([I]). Parameters: r = q = 0.01, uiq = 6, m = 5. 

225 of infected arrivals are shown in figures E] and HJ Figure [3] shows that the heterogeneous 
approximation (dashed lines) is in exact agreement with the agent-based model for the 
whole range of parameters p and w. The precision of the homogeneous approximation is 
also high for a large range of parameter values. As illustrated in figure HJ the absolute error in 
estimation of the disease prevalence of the approximation is maximal for intermediate values 
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230 of infectiousness p, but still less than 0.05. The only qualitative discrepancy between the 
approximation and the agent-based model emerges at low infectiousness p for zero infected 
arrivals (w = 0). Here, the homogeneous model predicts the existence of an epidemic 
threshold, whereas in the agent-based simulation and the heterogeneous approximation the 
disease is found to persist for any finite infectiousness p. 



235 V. ABSENCE OF EPIDEMIC THRESHOLDS 

Summarizing the results shown so far, we can say that ongoing epidemic dynamics quickly 
leads to the formation of networks with finite variance. Generally, one would expect that 
such networks should exhibit a finite epidemic threshold. Nevertheless, the heterogeneous 
approximation and simulations indicate that the epidemic can persist in these networks for 
240 any finite positive value of the infectiousness. Let us therefore investigate the apparent ab- 
sence of the threshold in greater detail. In the following we confirm this absence by what 
is essentially a less naive version of the naive line of reasoning presented in the introduc- 
tion. Throughout the argument we will only consider the case where all arriving agents are 
susceptible, as the epidemic threshold is harder to define otherwise. 

245 The epidemic threshold is commonly defined as the minimal value of the infectiousness, 
below which a randomly picked node is susceptible with probability 1. In any finite network 
this implies that below the epidemic threshold each individual node is susceptible. By 
contrast in the thermodynamic limit there can still be a finite number of infected nodes as 
long as the density of such nodes in the network is zero, i.e. [I] = 0. In the following we 

250 denote the state below the epidemic threshold as the disease-free state, but recognize that 
there may be still a finite number of infected agents. 

Because the density of infected vanishes in the disease-free state, it is apparent that the 
degree distribution must be independent of the parameters p and r. This can also be verified 
in the heterogeneous approximation. Furthermore, in the thermodynamic limit, i.e. N — > oo, 
255 the degree distribution converges to the power law Pk ~ k~ 3 for k > m and thus also the 
degree variance of susceptibles erf diverges in the disease-free state. Considering any state 
with stationary disease prevalence [/]*, mean degree (k)*, susceptible mean degree and (ks)*, 
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FIG. 4. Performance of the homogeneous approximation in estimating the disease prevalence 
[/]* in the (p,w) parameter space. Top-left: individual-based simulations. Top-right: homoge- 
neous approximation. Bottom: absolute difference between individual-based simulations and the 
homogeneous approximation. Homogeneous approximation performs well for a large range of in- 
fectiousness p and infected arrival fractions w. Parameters: r = q = 0.01, rriQ = 6,m = 5. 

the heterogeneous approximation implies 

±(k 2 ) =g( m 2 + 2m-^)+ - r ((k 2 ) - (k 2 s )[S}*) . (5) 

Using (k 2 ) = (kg)[S] + (kj)[I] it follows that a 2 ->■ oo iff a 2 s -> oo. In addition, the 
260 steady-state solution of the coarse-grained heterogeneous approximation provides that l/og 
vanishes only when [J] =0. Thus, the disease prevalence can vanish only when the degree 
distribution is scale- free. 

The argument above provides a substantiation for our intuition that, in the disease-free 
state, preferential attachment will eventually lead to the formation of scale-free topologies, 
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265 where the degree variance of the susceptibles diverges. Considered on its own this argument 
does not preclude the existence of an epidemic threshold, below which the network is disease- 
free and scale-free. Although we know that the epidemic threshold vanishes in models of 
epidemic spreading on static scale-free networks, the same is not necessarily true in the 
adaptive network considered here. We therefore explicitly compute the epidemic threshold 

270 by checking the stability of the disease-free state against invasion of the epidemic. 

In a dynamical system, a given steady state is stable if all eigenvalues of the corresponding 
Jacobian matrix have negative real parts. The Jacobian matrix is obtained by a local 
linearisation of the system around the steady state under consideration. For the disease- 
free state, we can obtain the Jacobian J directly from the coarse-grained heterogeneous 

275 approximation (equation (jHJ)), which yields 



/ 



-q + pm — r —p —p + 4pm 
2rm —q — 2r 2r 

ri r 2 r 3 



(6) 



where 



Ti = p + <jg + qm + 2mr 



1 2 = O a 

2m & 



r 

2 

Q 



P 2 , 

(To + r . 

2m s 2 



(7) 



We can now compute the epidemic threshold p* as the critical threshold in p where the 
disease-free state loses its stability through a transcritical bifurcation. Since the Jacobian 
J depends only on the value of the susceptible degree variance cr| but not on its partial 

280 derivatives, we can treat the susceptible degree variance <r| as a parameter and obtain the 
epidemic threshold as its function by numerical continuation of the bifurcation point. For 
sufficiently large degree variance <r|, we find the scaling relation p* ~ (k)/a 2 s (figure ED, 
which recovers a famous result from epidemics on static networks. The epidemic threshold 
thus vanishes in the thermodynamic limit when <t§ — > oo. 

285 The two paragraphs above provide a clearer picture of the paradoxical situation sketched 
already in the introduction. If the disease prevalence vanishes then the network becomes 
scale-free, however the scale-free network has a vanishing epidemic threshold and thus sup- 
ports finite prevalence. In principle this paradox can be resolved in two ways, a) either 
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FIG. 5. The epidemic threshold p* as a function of a given value of degree variance a 2 . Results 
have been obtained by numerical continuation of a bifurcation line from equation ([3]) . The epidemic 
threshold scales in the same form as on static networks, i.e. p* ~ (k)/a 2 , showing that the epidemic 
threshold vanishes in the thermodynamic limit. Parameters: r = q = 0.01, w = 0, mo = 6, m = 5. 

the network does not become stationary but switches back and forth between scale-free and 
290 exponential phases, or b) the network approaches a state where the epidemic persists in an 
exponential topology. 

It is apparent that in finite networks the temporal dynamics can lead to a disease-free 
scale-free state. In this case the network goes through an initial exponential phase where 
an epidemic threshold exists that leads to the extinction of the epidemic. Subsequently, 
295 a scale-free topology is build up, in which the epidemic threshold vanishes. However the 
epidemic cannot reappear as no infected are left which could reignite the epidemic. 

We emphasize that the mechanism above is a finite-size effect. In the thermodynamic 
limit there would generally be some infected left even in the disease-free state. To see 
this consider that infected are removed exponentially, so that after any finite time there 
300 will still be infected left which can launch the network back into an endemic state, such 
that scenarios a) and b) are the only possible outcomes. The finite-size extinction of the 
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epidemic is demonstrated explicitly in figure El which shows that the fraction of agent-based 
simulation runs in which the epidemic persists increases with the initial network size Mq. 
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FIG. 6. The fraction of surviving runs for low infectiousness p as a function of the initial network 
size Mq. The fraction of surviving runs increases as the initial network size mo increases. In order to 
ensure the same initial average degree m, the Barabasi- Albert growth model with m = 5 and tuq = 6 
is iterated until the network reaches size Mq. Then nodes are assigned states and the infection, 
removal, and network growth processes take place simultaneously. Parameters: r = q = 0.01, 
w = 0, mo = 6, m = 5. 



Precluding the finite-size effect described above, simulation runs in the parameter range 
305 considered so far approach a finite prevalence for any positive value of infectiousness. The 
paradox outlined above is thus resolved by alternative a) the model does not have an epi- 
demic threshold although the evolving topologies remain exponential. As a final step in 
our exploration we now search for evidence of the temporal resolution of the paradox, the 
alternative b). 

310 Up to now we have considered relatively small removal rates r. At sufficiently high 
removal rate, the evolution of the disease and the topology exhibits dynamics different from 
the convergence to a stationary state. Figure [7] shows a representative evolution at non-zero 
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FIG. 7. The dynamical evolution in the [/] — (&) plane for high removal rate r. At high infectiousness 
p, the observed behaviour resembles a homoclinic trajectory. The network grows from a healthy 
initial state until an infected eventually arrives and the disease quickly spreads over the network 
and the infected individuals subsequently die. Homogeneous approximation (red circles) fails to 
capture this behaviour and predicts a stable equilibrium. Parameters: q = 0.01, r = 0.10, niQ = 6, 
m = 5, w = 0.1. 



fraction of infected arrivals w. At the combination of high infectiousness p and high removal 
rate r, the dynamics resembles a homoclinic trajectory in the (k) — [J] plane. The disease 
315 spreads quickly over the network and it covers the whole population immediately. Then 
disease-induced removals dominate and the population becomes extinct until new healthy 
individuals arrive and the disease spreading restarts, which leads to cycles of population 
growth and collapse. 

When all arrivals are susceptible, i.e. w — 0, the epidemics in the finite population 
320 disappears entirely and the collapse-and-growth cycle cannot be completed. A disease-free 
scale-free network then emerges. Based on the arguments above one can suspect that the 
same behaviour cannot occur in the thermodynamic limit. Instead it is likely that the 
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observed dynamics forms part of a homoclinic cycle, where long phases of very low disease 
prevalence are disrupted by sharp outbreaks. Similar dynamics have for instance been 

325 studied in an adaptive- network model of cooperation among agents {54J]. In the present 
model, the moment closure approximation fails to capture this behaviour and predicts a 
stable steady-state solution for high values of r. 

The observed dynamics at high removal rates would correspond to diseases with very 
high mortality where infected individuals die almost immediately. This is reminiscent of the 

330 massive pandemics in history, where humanity was exposed to new pathogens with very high 
virulence. While such diseases could lead to the deaths of large fractions of populations, 
continuous supply of healthy individuals through immigration provided new hosts to the 
disease and introduced new bursts of disease spreading which caused repeated epidemics 
cycles. 



335 VI. DISCUSSION 

In the present paper we have investigated the dynamics of a fatal disease in a growing 
population. Our main finding is that no epidemic threshold exists in this model although 
the variance of the degree distribution remains finite. Thus even in evolving exponential 
network topologies "unlikely" diseases with very low infectiousness can persist indefinitely. 

340 We presented a detailed analytical exploration for the case of low removal rate, where the 
prevalence of the disease reaches a stationary level. In the growing population the ongoing 
epidemic dynamics eliminates the nodes of high degree and thus leads to the formation of 
exponential topologies for which the variance of the degree distribution is finite. However, 
this mechanism only reduces the width of the degree distribution so far that the epidemic can 

345 still persist. Thus for any finite value of the infectiousness the network adapts its topology 
such that the variance of the degree distribution is lower to point where the epidemic can 
still be sustained, which explains the observed absence of the epidemic threshold. 

When the removal rate is sufficiently high, we observed a much more dynamic picture, 
where degree variance acts as a supply of fuel that slowly builds up before being burned 

350 in large epidemic outbreaks. These dynamics are reminiscent of the massive deaths from 
epidemic diseases in newly forming urban areas. A detailed analysis thus appears as a 
promising target for future work. 
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The dynamical feedback between the population structure and the epidemic disease has 



been so far studied in a number of articles in the past five years |30|, l33M37j . However, models 
355 captured mostly social interactions in non-fatal diseases. A significant obstacle to progress 
in this line of work is that the network evolution in most models is driven by behavioural 
changes of individuals (whom to meet, how often to wash hands). However, for behaviour 
often no records exist such that model predictions cannot easily be compared with real 
world data. By contrast, in the model proposed here, the social network evolves due to 
360 demographic processes such as migration and death on which data may be easier to obtain. 

We believe that the proposed model is relevant for epidemics in rapidly growing cities, 
especially in the developing countries. We hope that in this context connecting the model 
to real world data will be feasible in the future. One can easily imagine extensions of 
the present model that incorporate policy measures such as vaccination, quarantine, or 
365 regulation of migration. We hope that this will in the future lead to the formulation of more 
efficient policies for combating epidemic diseases. 



[1] Hays, J.N. 2005 Epidemics And Pandemics: Their Impacts on Human History. Santa Barbara, 
California: ABC-CLIO. 

[2] World Malaria Report 2011, Geneva, Switzerland: WHO Press, World Health Organization. 
370 [3] Keeling, M. J. & Eames, K. T. D. 2005 Networks and epidemic models. J. R. Soc. Interface 
2, 295-307. 

[4] Bansal, S., Grenfell, B. T. & Meyers, L. A. 2007 When individual behaviour matters: homo- 
geneous and network models in epidemiology. J. R. Soc. Interface 4, 879-891. 

[5] Funk, S., Salathe, M. & Jansen, V. A. A. 2010 Modelling the influence of human behaviour 
375 on the spread of infectious diseases: a review. J. R. Soc. Interface 7, 1247-1256. 

[6] Moore, C. & Newman, M. E. J. 2000 Epidemics and percolation in small-world networks. 
Phys. Rev. E 61, 5678-5682. 

[7] Pastor- Satorr as, R. & Vespignani, A. 2001 Epidemic dynamics and endemic states in complex 
networks. Phys. Rev. E 63, 066117. 
380 [8] Moreno, Y., Pastor- Satorr as, R. & Vespignani, A. 2002 Epidemic outbreaks in complex het- 
erogeneous networks. Eur. Phys. J. B 26, 521-529. 



20 



[9] Pastor- Satorr as, R. & Vespignani, A. 2001 Epidemic spreading in scale-free networks. Phys. 
Rev. Lett. 86, 3200-3203. 
[10] May, R. M. & Lloyd, A. L. 2001 Infection dynamics on scale-free networks. Phys. Rev. E 64, 
385 066112. 

[11] Castellano, C. & Pastor- Satorr as, R. 2010 Thresholds for epidemic spreading in networks. 
Phys. Rev. Lett. 105, 218701. 

[12] Barthelemy, M., Barrat, A., Pastor-Satorras, R. &: Vespignani, A. 2004 Velocity and hierar- 
chical spread of epidemic outbreaks in scale-free networks. Phys. Rev. Lett. 92, 178701. 
390 [13] Pastor-Satorras, R. & Vespignani, A. 2002 Immunization of complex networks. Phys. Rev. E 
65, 036104. 

[14] Zanette, D. H. & Kuperman, M. 2002 Effects of immunization in small-world epidemics. 
Physica A 309, 445-452. 

[15] Cohen, R., Havlin, S. & Ben-avraham, D. 2003 Efficient Immunization Strategies for Computer 
395 Networks and Populations. Phys. Rev. Lett. 91, 247901. 

[16] Holme, P. 2004 Efficient local strategies for vaccination and network attack. Europhys. Lett. 
68, 908. 

[17] Meyers, L. A., Pourbohloul, B., Newman, M. E. J., Skowronski, D. M. & Brunham, R. C. 
2005 Network theory and SARS: predicting outbreak diversity. Journal of Theoretical Biology 
4oo 232, 71-81. 

[18] Newman, M. E. J., Forrest, S. &; Balthrop, J. 2002 Email networks and the spread of computer 
viruses. Phys. Rev. E 66, 035101 (R). 

[19] Hufnagel, L., Brockmann, D. & Geisel, T. 2004 Forecast and control of epidemics in a global- 
ized world. Proc. Natl. Acad. Sci. USA 101, 15124-15129. 
405 [20] Colizza, V., Barrat, A., Barthelemy, M. & Vespignani, A. 2006 The role of the airline trans- 
portation network in the prediction and predictability of global epidemics. Proc. Natl. Acad. 
Sci. USA 103, 2015-2020. 

[21] Read, J. M., Eames, K. T. D. & Edmunds, W. J. 2008 Dynamic social networks and the 
implications for the spread of infectious disease. J. R. Soc. Interface 5, 10011007. 
4io [22] Belik, V., Geisel, T. & Brockmann, D. 2011 Natural Human Mobility Patterns and Spatial 
Spread of Infectious Diseases. Phys. Rev. X 1, 011001. 

[23] Pastor-Satorras, R. & Vespignani, A. 2002 Epidemic dynamics in finite size scale-free networks. 



21 



Phys. Rev. E 65, 035 108 (R). 
[24] Eguiluz, V. M. & Klemm, K. 2002 Epidemic threshold in structured scale-free networks. Phys. 
4i 5 Rev. Lett 89, 108701. 

[25] Read, J.M. & Keeling, M.J. 2003 Disease evolution on networks: the role of contact structure. 

Proc. R. Soc. B 270, 699-708. 
[26] Serrano, M. A. & Boguna, M. 2006 Percolation and epidemic thresholds in clustered networks. 

Phys. Rev. Lett 97, 088701. 
420 [27] Mossa, S., Barthelemy, M., Stanley, H. E. & Amaral, L. A. N. 2002 Truncation of power law 

behavior in scale-free network models due to information filtering. Phys. Rev. Lett 88, 138701. 
[28] Grabowski, A. & Kosinski, R. A. 2004 Epidemic spreading in a hierarchical social network. 

Phys. Rev. E 70, 031908. 
[29] Newman, M. E. J. 2002 Spread of epidemic disease on networks. Phys. Rev. E 66, 016128. 
425 [30] Gross, T., D'Lima, C.J.D. & Blasius, B. 2006 Epidemic dynamics on an adaptive network. 

Phys. Rev. Lett. 96, 208701. 
[31] Gross, T. & Blasius, B. 2008 Adaptive coevolutionary networks: a review. J. R. Soc. Interface 

5, 259-271. 

[32] Gross, T. & Sayama, H. (Eds.) 2009 Adaptive Networks. Theory, Models and Applications. 
430 New York: Springer Verlag. 

[33] Shaw, L.B. & Schwartz, I. B. 2008 Fluctuating epidemics on adaptive networks. Phys. Rev. E 
77, 066101. 

[34] Risau-Gusman, S. & Zanette, D. H. 2009 Contact switching as a control strategy for epidemic 
outbreaks. Journal of Theoretical Biology 257, 52-60. 
435 [35] Graser, O., Hui, P.M. &; Xu, C. 2011 Separatrices between healthy and endemic states in an 
adaptive epidemic model. Physica A 390, 906-913. 

[36] Wang, B., Cao, L., Suzuki, H. & Aihara, K. 2011 Epidemic spread in adaptive networks with 
multitype agents. J. Phys. A 44, 035101. 

[37] Marceau, V., Noel, P.-A., Hebert-Dufresne, L., Allard, A. & Dube, L.J. 2010 Adaptive net- 
440 works: coevolution of disease and topology. Phys. Rev. E 82, 036116. 

[38] Shaw, L. B. & Schwartz, I. B. 2010 Enhanced vaccine control of epidemics in adaptive net- 
works. Phys. Rev. E 81, 046120. 

[39] Lagorio, C. , Dickison, M., Vazquez, F., Braunstein, L.A., Macri, P. A., Migueles, M.V., Havlin, 



22 



S. h Stanley, H.E. 2011 Quarantine-generated phase transition in epidemic spreading. Phys. 
44 5 Rev. E 83, 026102. 

[40] Kamenev, A. & Meerson, B. 2008 Extinction of an infectious disease: a large fluctuation in a 

nonequilibrium system. Phys. Rev. E 77, 061107. 
[41] Schwartz, I. B., Billings, L., Dykman, M. &: Landsman, A. 2009 Predicting extinction rates 

in stochastic epidemic models. J. Stat. Mech., P01005. 
450 [42] Hayashi, Y., Minoura, M. & Matsukubo, J. 2004 Oscillatory epidemic prevalence in growing 

scale- free networks. Phys. Rev. E 69, 016112. 
[43] Poncela, J., Gomez-Gardehes, J., Traulsen, A. &; Moreno, Y. 2009 Evolutionary game dynam- 
ics in a growing structured population . New J. Phys. 11, 083031. 
[44] Anderson, R. M. & May, R. M. 1991 Infectious Diseases of Humans: Dynamics and Control. 
455 New York: Oxford University Press. 

[45] Guerra, B. & Gomez-Gardenes, J. 2010 Annealed and mean-field formulations of disease 

dynamics on static and adaptive networks. Phys. Rev. E 82, 035101 (R). 
[46] Barabasi, A. & Albert, R. 1999 Emergence of Scaling in Random Networks. Science 286, 

509-512. 

460 [47] Keeling, M. J. 1999 The effects of local spatial structure on epidemiological invasions. Proc. 
R. Soc. B 266, 859-867. 

[48] Bauch, C. T. 2005 The spread of infectious diseases in spatially structured populations: an 

invasory pair approximation. Mathematical Biosciences 198, 217-237. 
[49] House, T. & Keeling, M. K. 2011 Insights from unifying modern approximations to infections 
465 on networks. J. R. Soc. Interface 8, 67-73. 

[50] Demirel, G., Vazquez, F. , B6hme, G. A. & Gross, T. Moment closure approximations of the 

adaptive voter model, In preperation. 
[51] Noel, P.-A., Davoudi, B., Brunham, R. C., Dube, L. J. & Pourbohloul, B. 2009 Time evolution 

of epidemic disease on finite and infinite networks. Phys. Rev. E 79, 026101. 
470 [52] Lindquist, J., Ma, J., van den Driessche, P. & Willeboordse, F. H. 2011 Effective degree 

network disease models. J. Math. Biol. 62, 143-164. 
[53] Gleeson, J. P. 2011 High-accuracy approximation of binary-state dynamics on networks. Phys. 

Rev. Lett. 107, 068701. 

[54] Zschaler, G., Traulsen, A. &; Gross, T. 2010 A homoclinic route to asymptotic full cooperation 



23 



475 in adaptive networks and its failure. New J. Phys. 12, 093015. 



